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We measured the local density of states of a quasi two-dimensional electron system (2DES) near 
defects, artificially created by Ar-ion sputtering, on surfaces of highly oriented pyrolytic graphite 
(HOPG) with scanning tunneling spectroscopy (STS) in high magnetic fields. At valley energies of 
the Landau level spectrum, we found two typical localized distributions of the 2DES depending on 
the defects. These are new types of distributions which are not observed in the previous STS work 
at the HOPG surface near a point defect [Y. Niimi et al, Phys. Rev. Lett. 97, 236804 (2006).]. 
With increasing energy, we observed gradual transformation from the localized distributions to the 
extended ones as expected for the integer quantum Hall state. We show that the defect potential 
depth is responsible for the two localized distributions from comparison with theoretical calculations. 

PACS numbers: 73.43.-f, 71.70.Di, 68.37.Ef, 71.20.Tx 



Two-dimensional electron systems (2DESs) exhibit fas- 
cinating quantum phenomena at low temperature. The 
quantum Hall (QH) effect in high magnetic fields and the 
Anderson localization are two well-known examples 
These phenomena have been studied mainly by means of 
transport measurements. However, it is generally diffi- 
cult to investigate the electronic states of the 2DESs in 
real space on nanometer scale with local probes, since 
they are usually formed at hetero junctions a few hun- 
dreds nanometers below the semiconductor surfaces. 

Recently, 2DESs formed at semiconductor surfaces 0, 
0, 0, [H and semimetal ones 0, 0] have been studied 
with scanning tunneling microscopy and spectroscopy 
(STM/STS). These are powerful techniques to investi- 
gate the local density of states (LDOS) at sample sur- 
faces. Morgenstern et al. [3[ observed clear Landau quan- 
tization as well as complicated patterns of the LDOS 
depending on bias voltage for a two-dimensional elec- 
tron gas (2DEG) at a cleaved InAs(llO) surface with 
submonolayer iron deposition. However, the localized 
and extended LDOS distributions are not clearly distin- 
guished in their measurements. 

More distinct localization and extension of the LDOS 
depending on energy, which indicate the possible QH 
state, have been observed near a point defect at sur- 
faces of highly oriented pyrolytic graphite (HOPG) Q- 
At valley energies of the Landau levels (LLs), a circu- 
lar distribution of the LDOS was observed near the de- 
fect. The distribution with a radius comparable to the 
magnetic length Ib (= \/h/eB where B is magnetic 
field) was semiquantitatively explained by the calculated 
LDOS for 2DEG in magnetic fields in the 1/r poten- 
tial In transport measurements for HOPG @, a 
Hall resistance plateau was observed, which is also in- 
dicative of the QH effect in this material. Unambiguous 
QH plateaus were observed in the 2DESs at a single layer 
of graphite (graphene) [10, LLl| and bilayer graphene [To[ . 
Since these thin-layer graphite systems including HOPG 



have the surface 2DEGs, they provide ideal arenas for in- 
vestigation of quantum phenomena in 2DESs, especially 
the QH state, with the STM/STS techniques. 

In this Letter, we studied the LDOS at HOPG surfaces 
with randomly distributed defects which are artificially 
created by Ar-ion sputtering. Two types of the localized 
LDOS distributions are observed at the valley energies of 
the LLs in the differential tunnel conductance (d//dV) 
images. As we increase energy, the distributions are ex- 
tended to follow the complicated potential landscapes. 
We reveal that the functional forms of the defect poten- 
tials are closely related to the LDOS distributions. 

The STM/STS measurements were performed at tem- 
peratures below 30 mK and in magnetic fields up to 6 T 
using an ultra low temperature STM [12[ . Electrochemi- 
cally etched W wires were used as STM tips. The dl /dV 
curves and images were taken by the lock-in technique 
with a bias modulation V mo a of 0.5 or 1.0 mV at a fre- 
quency of 412 Hz. The HOPG sample [13[ was cleaved 
in air and then quickly loaded into an ultra high vacuum 
chamber of the STM. Defects were made by sputtering 
the sample surface in situ with a 30 eV Ar-ion beam for 
a few seconds in an Ar atmosphere of 10 -6 Pa. Higher 
energy beams for a longer time made the surfaces too 
rough for the STM/STS measurements. 

Figure 1(a) shows a typical STM image over 100 x 100 
nm 2 at the sputtered HOPG surface. 16 defects (De- 
fects 1 ~ 16) are seen in this image. The averaged defect 
density is 2 x 10 11 cm -2 . The defects have nearly the 
same diameters (6^8 nm) which are much larger than 
that of the point defect (< 1 nm) studied in Ref. 0- 
They are randomly distributed on the surface. In Fig. 
1(b), we show close-up images of two of them, where 
complicated inside structures are clearly seen. Just out- 
side the defects, the (V^ x y/S)R30° superstructure and 
honeycomb one consisting of the B-site carbon atoms are 
observed [3 El. 

In Fig. 1(c), we show tunnel spectra at several different 
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FIG. 1: (a) Typical STM image at the sputtered HOPG sur- 
face [100 x 100 nm 2 , V = 120 mV, I = 0.2 nA]. The numbers 
are assigned for all (sixteen) defects in this image, (b) Close- 
up STM images [10 x 10 nm 2 ] of Defects 1 and 4. The diamond 
and honeycomb represent the (y/3 x \/3)R30° superstructure 
and honeycomb one. (c) Top figure shows tunnel spectra at 
the sputtered surface in several different magnetic fields per- 
pendicular to the graphite basal plane. Each spectrum is 
shifted for clarity. Bottom figure shows the LL peak energies 
as a function of B. The dotted and broken lines are linear and 
square root fittings of the experimental data, respectively. 



magnetic fields averaged over the whole scan area of Fig. 
1(a). Clear LL peaks similar to those in the previous 
works jf| 0] are observed. The Landau index (n) for 
each peak is assigned as in Ref. Q . By comparing with 
the calculated surface LDOS, the effective thickness of 
this area is estimated as about 20 layers (= 6.7 nm) of 
graphene [6]. Note that the Fermi energy (Ep) here is 
shifted by +30 meV, compared to those in the previous 
works 0, 0] • This is because a large number of defects 
made by sputtering inclined the electron and hole balance 
to excess electrons. 

It is intriguing to recognize that the LLs have two dif- 
ferent magnetic field variations depending on n. For ex- 
ample, the peak energies of LL1, LL2 and LL3 near Ep 
are proportional to B as expected for the conventional 
2DEG. Meanwhile, the LLs far from Ep denoted by open 
symbols in Fig. 1(c) have a \[B dependence. If we fit the 
latter LLs with E n = 2ehv 2 F B\n\, the formula of the LL 
energies for massless Dirac Fermions [16j . we obtain the 
Fermi velocity vp of 1 x 10 6 m/s. This value is consistent 
with the estimation from the transport measurements of 



graphene 1(J HI- These facts, which are consistent with 
similar results of Ref. [17[ , indicate that there are linear 
and parabolic subbands at the HOPG surface 0. 

Next, in Figs. 2(a)-2(n), we show dl/dV images at 
different bias voltages at 6 T in the same scan area as 
Fig. 1(a). Note that the same contrast is used for all 
the images. At the valley energy in between LL0, —1 
and LL1 [Fig. 2(a)], the electronic states are strongly or 
weakly localized near some of the defects. The localiza- 
tion is gradually dispersed with increasing energy until 
LL1 [Figs. 2(b)-2(i)]. 

First of all, we focus on the dl/dV image at the valley 
energy in Fig. 2(a). We found two new types of localized 
LDOS distributions depending on the defects, which are 
not observed in the previous STS work 0. One is a 
distribution where the LDOS amplitude has a maximum 
just on defect (Type I: denoted by cross). The other 
is a ring distribution around defect (Type II: denoted 
by solid circle). We also found that five defects do not 
support any localized states nearby (denoted by triangle). 
In Figs. 2(o) and 2(p), we show cross sections of Type I 
(Defect 1) and Type II (Defect 5) at different fields. Both 
distributions have diameters of the order of Ib (oc 1/y/B) 
at each field, but the peak amplitude of Type I is much 
larger than that of Type II. The latter is qualitatively 
different from the LDOS distribution observed near the 
point defect Q, while the former is seemingly similar to 
it. However, Type I does not have an appreciable satellite 
ring unlike the point defect case. And it suddenly changes 
to Type Il-like in between 2 and 3 T with decreasing field. 
Such behavior is also observed for Defects 2, 6, and 16 
but not for Defects 8 and 9. 

As we increase energy [Figs. 2(b)-2(i)], the LDOS lo- 
calization is gradually dispersed, showing complicated 
extended patterns. It should be noted that the electronic 
states are extended to keep away, for instance, regions 
denoted by arrows (Regions A and B in Fig. 2(h)). In 
contrast, as energy exceeds LL1 by 2 meV [Fig. 2(k)], 
the electronic states are existing in those regions, and 
also Type II distributions can be seen near some of the 
defects (center- left). A similar evolution is observed at 
energies in between other two successive LLs below LL3 
[Figs. 2(m) and 2(n)]. 

The localized and extended states can be discriminated 
more quantitatively by plotting the coefficient of varia- 
tion (A) for each dl/dV image (which is given by dividing 
the standard deviation by the average) as a function of 
bias voltage in Fig. 2(q). It is obvious from this plot that 
the spatial variations of LDOS are larger at the valley en- 
ergies in between two LLs, while they are smaller at the 
peak energies of the LLs. This is an additional indication 
for a transition from localized to extended states. 

Figure 2(r) is the dl/dV image at the peak energy 
of LL0, —1. As is shown in the previous STS works on 
graphite [6, 7], this peak is much more sensitive to the 
local electrostatic potential compared to the other LL 
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FIG. 2: (a)-(n) dl/dV images in the same scan area as Fig. 1(a) at various bias voltages [100 x 100 nm 2 , B = 6 T, T — 30 
mK, / = 0.2 nA, V mod = 1.0 mV]; (a) -0.1, (b) 1.1, (c) 2.3, (d) 3.5, (e) 4.7, (f) 5.9, (g) 7.1, (h) 8.3, (i) 9.5, (j) 10.7, (k) 11.9, 
(1) 13.1, (m) 14.3, (n) 31.6 mV. The defect positions are denoted by the crosses, solid circles and triangles, depending on the 
kind of the localized LDOS distribution (see text). The arrows in (h) and (k) show the regions (Regions A and B) where the 
electronic states are hardly extended with increasing energy, (o), (p) Cross sections of the dl/dV images averaged radially near 
(o) Defect 1 and (p) Defect 5 at the valley energies in between LL0, — 1 and LL1 at different fields. Each cross section is shifted 
by (o) 2 or (p) 0.4 nA/V for clarity. The arrows in (p) show the approximate diameters of the LDOS rings, (q) Average (solid 
circle) and coefficient of variation (open circle) of dl/dV in the whole scan area as a function of bias voltage, (r) dl/dV image 
at the peak energy of LL0, — 1 ( — 12 mV). (s) Tunnel spectra averaged over the whole scan area (black), over 10 x 10 nm 2 near 
Defect 1 (blue), Defect 5 (red), Region A (green), and Region B (pink). Each spectrum is shifted by 1 nA/V for clarity, (t), 
(u) Schematics of the wave functions for the ground states in (t) the 1/r and (u) harmonic potentials. 



peaks. This is attributable to much narrower but dis- 
persive subband structure of LL0, —1 along the k z axis 
[see Fig. 2(b) of Ref [6|]]. Thus, we assume that the im- 
age at the LL0, —1 peak energy is indicative of the sur- 
face potential distribution which should be crucial to the 
electron localization. For instance, there are not strong 
potential variations around the defects which do not sup- 
port any localized states in Fig. 2 (a). Figure 2 (s) shows 
the tunnel spectra near Defect 1 (Type I) and Defect 
5 (Type II). The LL0, — 1 peak energies near Defect 1 
(V = -18 mV) and Defect 5 (-15 mV) are shifted to 
the negative bias voltage side compared to that averaged 
over the whole scan area (—12 mV). This indicates that 
the potential of Defect 1 is deeper than Defect 5 and that 
both are confining ones. The same tendencies are seen 
for the other defects for Types I and II. On the other 
hand, in Regions A and B, the LL0, —1 peak energies are 
shifted to the positive energy side [Fig. 2(s)], which in- 



dicates that the potentials of those regions are relatively 
high compared to the other regions. These results sug- 
gest that the observed energy and spatial variations of 
the LDOS are determined by the surface potential distri- 
bution. This is expected behavior for the localized and 
extended states in the QH effect [19]. 

The two different localized LDOS distributions, i.e., 
Types I and II, can be qualitatively explained by recent 
calculations on eigenfunctions of simple 2DEG in mag- 
netic fields [8]. With a spatial variation of electrostatic 
potential, the degeneracy of the LLs is lifted. In the 
case of 1/r potential, the angular momentum (L z = £ z h) 
of the ground state is zero for each LL. Thus, its wave 
function has a maximum at the origin to minimize the 
potential energy [Fig. 2(t)]. On the other hand, the 
ground state in the harmonic potential has t z = n, and 
the maximum of electron probability forms a ring around 
the origin [Fig. 2(u)]. The distributions of Types I and II 
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FIG. 3: (a)-(g) dl/dV images near a defect at various bias 
voltages [50 x 50 nm 2 , B = 6 T, T = 30 mK, I = 0.2 nA, 
V mod = 0.5 mV]; (a) 29, (b) 30, (c) 32, (d) 33, (e) 34, (f) 
35, (g) 37 mV. The same contrast is used for all the dl/dV 
images. The inset of figure (c) shows the same dl/dV image in 
color scale, (h) Cross sections of the dl/dV images averaged 
radially near the defect. The inset is calculated LDOS at 
different energies in the 1/y/r 2 + d 2 potential shown in Fig. 
11 of Ref. [8]. 



have thorough theoretical explanations for this observa- 
tion at the moment. A qualitatively similar energy de- 
pendence of localized state has been given in the previous 
calculations with the 1/^/r 2 + d 2 potential [inset of Fig. 
3(h)]. The essential point here is the spectral weight con- 
version from the ground state to the higher-order excited 
states. The improvement of the calculations by tuning 
the functional form of potential is highly desirable. 

In conclusion, we have measured the LDOS near de- 
fects artificially created by Ar-ion sputtering at HOPG 
surfaces in high magnetic fields with STM/STS. At the 
valley energies in between the adjacent LLs, the two new 
localized distributions were observed depending on the 
kinds of the defect potentials. With increasing energy, 
the localized distributions are extended to follow the po- 
tential landscapes created by the randomly distributed 
defects as is expected for the QH system. 

We thank D. Yoshioka for valuable discussions. This 
work was financially supported by Grant-in- Aid for Sci- 
entific Research on Priority Areas (Grant No. 17071002) 
from MEXT and ERATO Project of JST. Y.N. acknowl- 
edges the JSPS Research program for Young Scientists. 

note added: After submitting, we became aware of a 
related STM/STS work on a semiconductor surface [20| 
showing similar results to ours. 



are basically consistent with the ground state ones calcu- 
lated for the 1/r and harmonic potentials, respectively. 

There are, however, a few experimental observations 
which cannot be described well by theoretical calcula- 
tions. For example, the LDOS satellite ring for Type I 
is less obvious than that calculated for the 1/r potential. 
This can be alleviated with shallower potentials such as 
1/y/r 2 + d 2 (d is the depth from the surface). The local- 
ized LDOS distribution changes from Type I to II below 
3 T, while the calculated LDOS has Type I like distribu- 
tion down to much lower field [8[. 

Finally, let us describe a detailed energy dependence of 
Type I localized distribution in Fig. 3. This observation 
became possible near a similarly created but isolated de- 
fect at an HOPG surface with a less defect density. At 
the valley energy in between LL1 and LL2, the electronic 
state is localized just on the defect. With increasing en- 
ergy, however, the LDOS maximum is shifted from the 
center to the outside, and forms the ring- type distribu- 
tion [Figs. 3(a)-3(c)]. The LDOS ring gradually expands 
and disappears at the peak energy of LL2 [Figs. 3(d)- 
3(g)]. This evolution is sumarized in Fig. 3(h) where 
cross sections of the images are plotted. Similar evolu- 
tions are also but less obviously observed in the LDOS 
near other defects shown in Figs. 2(a)-2(i). We don't 
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